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ABSTRACT: 

It is shown that ion collisions introduce a slight reduction in the 
plasma frequency, along with an exponential decay of transient 
electron oscillations. The critical frequencies for penetration of a 
homogeneous ionosphere, for both isotropic and anisotropic iono- 
spheres, are determined. The characteristic waves , for electromagnetic 
propagation within a homogeneous anisotropic ionosphere, are developed 
by considering an infinite series of electron velocities, produced by 
Ampere's force law, reacting with the electron oscillation produced by 
an exciting electric field. The complex indices of refraction are 
determined, both from a dispersion equation and from a derivation of 
the Appleton equation, which uses an arbitrary selection of the co- 
ordinate axes, thus emphasizing the invariance of the Appleton equation. 
Vector and tensor algebra is used throughout the analytical developments. 
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Symbols 



E 0 = Vector magnitude of incident field 

B 0 = Geodesic field vector 
)c 0 = 6 o/c = Free space wave number 
' = Relative permeability constant 

€ '= Complex relative dielectric constant 

a 

€ ' = Complex relative dielectric tensor 



n = Complex vector index of refraction 

y = j )c e n = Complex vector propagation constant 
0= Angle between yand B 0 
5 = Cos • a 2 

X, y, Z = Cartesian coordinates 
r, 0,(p= Spherical coordinates 
£= Small displacement along X-axis 

F = Force vector 

e = Positive magnitude of electron charge 
m = Mass of electron 
- e/m = -1.77 x 10^ coul/Kg 

p = Average number of electron-ion collisions per second 
r = Radius vector 



V = Vector velocity of an electron 
N = Density of electrons per cubic meter 

_7 

ju 0 = 4 tt x 10 h/m = Permeability of free space 
1 ~9 

€ 0 « x 10 fd/m = Permitivity of free space 



w 



36tt 
N e 2 
m ^ 



= Resonant angular frequency, lossless plasma 



eB 0 = Cyclotron, or gyro, frequency 
m 

E, H = Electromagnetic field components 
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Symbols (Continued) 



D = Electric flux density vector 

i = Current density in amps, per sq. m. 
yO s = Surface charge density in coul. per sq. m. 
p = Charge density in coul. per cubic m. 

X = co p 2 / co 2 

Y - W b /co- J y 2 + Y 2 + Y 2 

x y 2 

Y = Y cos 0, Y = Y sin 0 cos <o, Y = Y = Y sin 0 sin <a= Y cos ib 

X V 2 L r 

~ ^y 2 + y 2 - ^ y 2 - Y 2 = Y s in 0 
X Y L 

r = i-j p/co 

M 2 = n 2 - 1 = n . n -1 

a . = Unit vectors 
i 

A . = a x a = a. sin r 
1 0 z 1 b 

*2 = ,a- e x l z ) l z = J 2 sin ? 

x = i - r x 




Y = X Y 




Z= 1 - x_ 

r 

. -1 E 

= tan y 

E 

X 

co = Critical frequency 
c 

co = Ordinary wave critical frequency 

co c ^ = Left circularly polarized wave critical frequency 

co D - Right circularly polarized wave critical frequency 

a = Lorentz conductivity tensor 

= Lorentz resistivity tensor 



E = V 

xy 



2 2 
E + E 

x y 



xy 
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Symbols (Continued) 



I = Identity tensor 

s = Length of path through the ionosphere 

g 

c « 3 x 10 m/sec = Speed of light in a vacuum 

V = Group speed 

g 

V = Phase speed 
P 

a = Attenuation constant 

= Phase constant 

0 O = Angle of incidence 

0^ = Angle of refraction 

R = Measure of polarization 

n.. = n, n. = Components of refractive tensor 
13 1 3 
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1. PLASMA FREQUENCY 
1. 1 Lossless Case 

Consider a lossless plasma having a density N (^f C ) 



! try 



0 f f ree electrons. If the electrons are compressed 
X laterally through a small displacement, £ , along the 
X-axis, a surface charge density 



X =0 X= i 
Figure 1 



-P s = -Ne i 



(1) 



may be considered to exist at y = £ . Then, from the definition of the 
electric displacement vector, from \= o to \= £ , 



D = 



Hence , 



D — f 0 E — Ne £ , 
X X 



or the force F is given by, 
X 



F = -eE = - Ne A 

X X ~ — 

^0 



Since force is given by the product of mass with acceleration , 

m S + Ne 2 £ = 0. 
dt 2 ~ 

2 2 

Substituting oo = Ne into (5), the harmonic equation for free oscill- 
P m 

ations of the electrons is obtained, hamely, 

d 2 £ 2 

-4- + “p «• 

dt^ F 



( 2 ) 

(3) 

(4) 

(5) 



( 6 ) 



The solution of (6) is sinusoidal with an angular frequency co . Thus oo 

P P 



is called the resonant angular frequency for the lossless plasma. 



1. 2 Lossy Case 

Suppose there are an average number, i> , of collisions per second 
between electrons and positive ions. This introduces an average loss of 
momentum, or an equivalent fluid resistance term, into equation (5), which 
equation becomes, 

2 2 
m d £ + pm d 4 + Ne £ = q. (7) 

dt 2 dt c 0 

Dividing by m, and using the operational method of substituting 



1 



into equation (7), the determinantal equation is obtained, 

(p 2 + vp + co 2 ) £ = 0. (8) 

r*' 

Equation (8) yields the characteristic roots, 

p = -v/2 ±j w p V i - ( WQ. (9) 

From equation (9), it may be concluded that ionic collisions in the 
plasma introduce a time decay constant of 2/v, and reduce the plasma 
frequency by a factor 4 ^ . For a low loss plasma, the 

latter factor usually is ignored. 



2. HOMOGENEOUS 
2 . 1 Low Loss Case 




electron is given by 



ISOTROPIC IONOSPHERE 

Consider a small collision frequency for a homo- 
geneous isotropic ionosphere of electron density 
N( e * ) . Let a uniform plane wave front, E 0 ge Ja,t , 
be incident upon the ionosphere at an angle of 
incidence 0 O . If r is a radius vector describing the 
electron having a velocity V, the force F upon an 



or 



F = -eE = 91 d^ r + v m _d r 

d, 2 dt 

d V + v V = -e_ E . 
dt m 



(10) 

(ID 



Substituting V e )Wt into (11) and solving for the steady state velocity, 



V = - e E 



= - e E 



( 12 ) 



m (jtd- y) jm( to -j v) 

The equivalent current density due to the motion of electrons then becomes. 



t = p V = -NeV = Ne 2 E 

m (v+ jto) 

Equations (12) and (13) show that the electron path, velocity, and current 
density are linear in the direction of E. 

In order to see the effects upon the electromagnetic wave. Maxwell's 
equations in the steady state phasor form will be used. 
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The Maxwell's equations are, div B = 0, div D = 0, curl E = - _3 b , 

dt 

curl H = 7 + 5D , B = ju' ju o H , 6 = f ' e 0 E . (14) 

Bt 

The corresponding complex form is. 



v. H = 0 , v .E = 0, v x E = -jco ji 0 H, 
v x H = l + j co E. 



(15) 



Substituting from equation (13), 



2 - 



vx H = Ne E + jco e 0 E = jcoe 0 # 
m {u + j co) 

[1 + ^ e2 — ] E , 



(16) 



or. 



j to e 0 m(p+ jo>) 
vxH = jco e 0 [ 1 + - 



co p 2 (v- jco) 
iuiv* + to 2 ) _ 



E = 



uo e c 



v co, 



( 1 ~ 0 “j 

ir + or oo( v + co z ) 



'P—T, 



(17) 



2 2 

in which co = Ne was also substituted in the right member of equation 
P m e 0 

(17). 

obtained , 



From equation (17), an equivalent complex dielectric constant c is 

c 



COr 



r '-j e" = ( 1- 
c 2 2 
CO +v 



-) -j 



_ i P CO 



2 2 
U + V L 



(18) 



Returning to euqations (15), the wave equation will be determined 
by elimination. 

- - 9 - 2 - 

Vx (vxE) =v(v. E) -V E = -V E = -jcOjj 0 A- 

V x H . 



(19) 



Using y = <*+ j j3 = jk 0 V c ' , and k 0 = co'. ju 0 c 0 = co/c , upon 

c 

substituting for vx H, equation (19) becomes, 

v 2 E + k 2 € ' E = 0. 

0 c 



(2 0) 



For v rsO , n ='J~P will be evaluated. 

n = U'-H') 1/2 ='^7 



// 1/2 

( 1-j _£_) 



V c* ( 1-j -i — ) = V e 
2 e' 



-j-* 



2,/V 



( 21 ) 
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